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To implement the BB84 decoy-state quantum key distribution (QKD) protocol over 
a lossy ground-satellite quantum uplink requires a source that has high repetition 
rate of short laser pulses, long term stability, and no phase correlations between 
pulses. We present a new type of telecom optical polarization and amplitude mod- 
ulator, based on a balanced Mach-Zehnder interferometer configuration, coupled to 
a polarization-preserving sum-frequency generation (SFG) optical setup, generating 
532 nm photons with modulated polarization and amplitude states. The weak coher- 
ent pulses produced by SFG meet the challenging requirements for long range QKD, 
featuring a high clock rate of 76 MHz, pico-second pulse width, phase randomization, 
and 98% polarization visibility for all states. Successful QKD has been demonstrated 
using this apparatus with full system stability up to 160 minutes and channel losses 
as high 57 dB 1 . We present the design and simulation of the hardware through the 
Mueller matrix and Stokes vector relations, together with an experimental implemen- 
tation working in the telecom wavelength band. We show the utility of the complete 
system by performing high loss QKD simulations, and confirm that our modulator 
fulfills the expected performance. 
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I. INTRODUCTION 



Utilizing ground-to-space quantum communications with satellites to achieve long- 
distance quantum key distribution (QKD) has been theoretically studied and experimentally 
proven feasible^. Ultimately, a global quantum network can be developed by using satel- 
lites as trusted nodes. An satellite uplink scheme is appealing because the energy-demanding 
and relatively complex photon source remains at the ground, easing the requirements on the 
satellite and granting the ability to use a variety of different sources. Practical QKD appa- 
ratuses use either weak coherent pulse (WCP) sourced^ED, or sources of entangled photon 
paird^HEO The major challenge of an uplink approach is the additional loss stemming from 
atmospheric turbulence, leading to low signal compared to the noise of detector dark counts 
and stray light. 

To perform QKD under such challenging conditions, it is desirable to have the source emit 
single photons within the shortest time window possible. This can be achieved with a mode- 
locked laser, which intrinsically provides short pulses and a high repetition rate. However 
such a laser has the inevitable problem of possessing phase correlations between consecutive 
pulses. This phase correlation violates an assumption of QKD security proofd^^, and is 
unsuitable as a photon source for QKD without sufficient precaution. To solve this problem, 
our source produces green photons at 532 nm wavelength by the sum-frequency generation 
(SFG) process, pumped by a pulsed 810 nm laser and a continuous wave (CW) laser at 
1550 nm. The hybrid nature of this design enables us to exploit the repetition rate and 
pulse width of the mode-locked Ti:Sapphire laser (810 nm), while phase randomization is 
accomplished with the short coherence length of the telecom laser. Importantly, the photons 
at 532 nm allow the use of detectors with the highest detection figure of meritP^ available. 
Additionally, we gain access to fast and stable modulation components designed for the 
telecom-band around 1550 nm. In this work, we mainly present the telecom intensity and 
polarization modulator (IPM) to encode quantum keys, as well as the QKD results to verify 
its performance as one of the building blocks. The SFG nonlinear optical setup is discussed 
elsewher#^. 

The desired performance for satellite-based QKD previously discussed in, e.g., RefsP 
ancPO, requires that we consider the following factors when implementing BB84 QKD 
protocol with the IPM: (1) Polarization encoding in the four non-orthogonal BB84 stated. 
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The main advantage of polarization encoding is that Earth's atmosphere preserves polar- 
ization with high fidelity^. (2) Compatible with the SFG setup to produce 532 nm green 
photons. (3) Compatible with Ti:Sapphire repetition rate, and readily upgradable to several 
GHz repetition rate. 

Previous implementations of polarization modulation for QKD include multiple laser 
diode^S, a single laser diode with multiple optical amplifier^, a single phase modulator with 
Faraday mirroi^, as well as a single polarization modulated. The first two approache d 10 1 25 1 , 
due to their using distinct sources of light, face the difficulty of making all output quantum 
states identical in frequency, bandwidth, and intensity; distinguishing information leads to 
security loopholes. The latter two method J^Sl will suffer polarization mode dispersion when 
two orthogonal modes have to propagate along unsymmetrical optical axes of the same 
modulator crystal or polarization maintaining fiber, leading to poor polarization states. 
Long-term thermal stability is also problematic for these implementations, and all exhibit 
pulse lengths at least two orders of magnitude longer than those obtained with a mode-locked 
laser. 

To avoid the above problems, we present our high speed polarization modulator in the 
balanced Mach-Zehnder Interferometer (MZI) configuration, including detailed descriptions 
of its design principles, and analysis of the experimental performance of the modulator. The 
paper is organized as follows: first we present an overview of the polarization and intensity 
modulator design, and its role in the entire QKD system; then we discuss the mathematical 
modelling of the telecom polarization and intensity modulator; finally, we show experimental 
results and the expected performance of decoy-state QKD based on our design. 

II. ALL-FIBER TELECOM BAND INTENSITY AND POLARIZATION 
MODULATION SYSTEM: DESIGN, MODELLING AND 
CHARACTERIZATION 

A. Overview of the QKD source design 

Implementing the BB84 protocol with decoy states requires that the output photons have 
at least two levels of average photon number. In the simplest case, this entails a signal state 
with average photon number /x and a decoy state with average photon number v. For both 
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levels, the modulator should output one of four polarization states (horizontal, vertical, 
diagonal, or anti-diagonal), chosen randomly. For all states, the average photon number per 
pulse is less than 1 to keep the multi-photon probability low 27 . However, the decoy-sate 
QKD protocol allows for a much larger average photon number (e.g. \i = 0.5) for signal 
pulses, as compared to WCP QKD protocols not employing decoy states^. 

Our QKD system is based on the SFG photon scheme as showed in Fig. [TJ in order to 
achieve the required modulation performance, pulse length, and phase randomization at the 
desired wavelength. The scheme includes a 1550 nm CW laser modulated in amplitude and 
polarization before being converted to 532 nm through SFG with strong 810 nm pulsed light 
from the mode-locked Ti:Sapphire laser. An overview of the SFG solution is illustrated in 
Fig.[j](a). The SFG green photons' field intensity 7 3 for Type-I PPKTP is proportional to the 
product Iil 2 where Ii and 7 2 represent the telecom and 810 nm optical power intensities. 
The linear dependence of the output intensity of green photons (whose spectrum is shown 
in Fig. [j] (b)) on the input intensities of both the telecom and 810 nm pumps allows direct 
modulation of the output light through modulation of only the telecom input. Since we 
are on at the regime of weak conversion, the spectrum of the generated photons will not be 
perturbed by the intensities of I\ and 7 2 - 

The QKD system has two basic parameters that determine the final key rate and max- 
imum distance: quantum bit error rate (QBER) and sifted key rate. A nonzero QBER 
is mainly caused by detector dark counts and polarization state visibility. To minimize 
QBER, by utilizing temporal filtering of detector and background noise, we incorporated 
thin-silicon SPAD detectors in Bob's receiver. The efficiency of these detectors peaks at a 
wavelength around 532 nm. 

The sifted key rate is determined by the repetition rate of the source, the average photon 
number /x, and channel losses. Pockels cells are capable of switching at such short wave- 
lengths, but they require driving electronics operating at a few kilovolts for a tt phase shift 
(14), preventing such modulators from being operated at the repetition rate of a Ti:Sapphire 
laser (^76 MHz), let alone approaching GHz rates desirable of high-speed QKD. Instead, our 
scheme in Fig. [j] allows for using conventional electro-optical (EO) waveguide modulators 
to match the desired repetition rate. 
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FIG. 1. (a) Generation of QKD polarization states at 532 nm by polarization-preserving SFG 
through a walk-off compensated PPKTP crystal pair (labeled as H and V, respectively). Light at 
810 nm produced by a Ti: Sapphire laser, mode-locked with 76 MHz repetition rate, is combined with 
1550 nm light with polarization state engineered by the IPM synchronized to the repetition rate of 
the TkSapphire laser. The SFG green (532 nm) photons combine properties from both pumping 
sources, including pulse shape and repetition rate from the Ti:Sapphire laser, and average photon 
number per pulse and polarization states from the modulated telecom laser, (b) The measured 
spectrum of green photons produced by our polarization-preserving SFG optical setup. Horizontally 
and vertically polarized photons have identical spectra. 



B. Overview of the telecom all-fiber modulation system 



The telecom modulation system for implementing decoy-state BB84 protocol consists 
of an intensity modulator followed by a polarization modulator, as realized by two phase 
modulators in the balanced MZI configuration that is illustrated in Fig. [2] (a). The physical 
apparatus of this scheme is pictured in Fig. [2] (b). The phase and intensity modulators are 
customized EO LiNbO x modulators from EOSpace, featuring low insertion loss (<2 dB) and 
external 50 Q termination of the radio frequency (RF) driving voltage. Their is as low 
as only a few volts, compatible with regular analog- digital mixed signal electronics. 

To ensure enough key is being received even when channel loss is high, the system should 
be capable of modulation at a few hundred MHz. The factors that limit the system clock 
rate are the finite bandwidth of the modulator drivers, the finite speed of field programmable 
gate array (FPGA) circuit board, and the interface circuits. The detailed design is displayed 
in Fig. 0(a). 
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The modulator driver circuit is comprised of digital-to-analog interface to bridge FPGA 
board and the power amplifier that are both adapted to be DC coupled to the RF ports of 
all modulators. To accomplish wide-band uniformity of the modulation, the driving signal 
bandwidth ranges from zero to several hundreds of MHz, resulting in reliable switching 
between any two polarization states. This is particularly important when the RF driving 
is in the lower frequency regime (e.g. during source testing), because the resulting heat 
accumulation is detrimental to the stability of polarization states, leading to poor QBER 
performance. 



C. Modelling of the telecom all-fiber modulation system 



We model the action of the modulation system using the Jones matrix approach assuming 
fully polarized input light and define the following Jones vector for the horizontally polarized 
state: 



H 







The system is comprised of an intensity modulator (IM), followed by a polarization beam 
splitter (PBS) with its fast axis rotated 45° from the input fast axis (represented by a 7r/4 
polarization rotation). Both input and output ports of the PBS are coupled to PM fiber. 
Following this is the balanced MZI, and after which is the polarization beam combiner (PBC) 
with PM input and single mode (SM) output. The Jones calculus is employed to integrate 
all of the optical components into a single matrix [J™ 

[J] = [QWP [PBC] [MZI] [PBS] [Ro Q + <&)] . (2) 

The last [QWP (— f )] is used to rotate from modulation in the circular-diagonal polarization 
plane to the linear polarization plane, and is performed by FPC_2 in the apparatus. 

We consider the intensity modulator separately as it acts globally on the total out- 
put intensity independent of polarization. Its mathematical model, as usually defined, is 
IM (Vb) = \ jl + &cos (J^tt + $1^ | where < b < 1 is modulation depth and $1 is the 
zero voltage phase, adjustable by applying an external DC bias. The applied voltage is Vo 
and the voltage parameter is defined in RefP^ as: V n 



— 4- in which r ? - 7 - is the linear 



electro-optical tensor element, no is the zero voltage refractive index, 7 is the optical confine- 
ment coefficient, L is the length of waveguide, and d is the gap distance of transmission line 
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FIG. 2. (a) Diagram of the intensity and polarization modulator. Signal and decoy states are 
realized by the intensity modulator, and polarization modulation by the balanced MZI, where each 
arm contains an EO phase modulator. The input and output to the modulator are a polarization 
beam splitter and a polarization beam combiner, respectively. The FPGA logic and digital Sz 
analog mixed signal board provide the encoding and driving voltages to the phase modulators 
and intensity modulator. The FPGA board is interfaced to three independent channels of fast 
digital to analog converter (D/A) circuits. The RF driving amplifier (DRV) are DC coupled to 
three ports of EO modulators. Two paddle-wheel fiber polarization controllers (FPCs) are used 
to adjust the input and output polarization states. Thick yellow lines (color) illustrate standard 
single mode fiber (SMF); thick blue lines (color) represent polarization-maintaining fiber (PMF); 
each red cross (color) indicates narrow key FC/PC connectors, enabling high-precision polarization 
coupling between two connecting optical fibers. The FPGA system clock is externally obtained 
from the PIN diode pulses of a mode-locked Ti-Sapphire laser. Those pulses are amplified via a 
pulse shaping circuit to generate 76 MHz TTL clock signal, (b) Photo of modulators including 
intensity (external PM fiber coupled) and phase modulators, as well as the polarization beam 
splitters and combiners (sealed in a thermally insulating box), (c) The core of this polarization 
modulator, comprised of a pair of phase modulators at the lower level, and one each of polarization 
beam splitter and combiner. 



that is excited by the modulator electrode input voltage. For the intensity modulator, the 
theoretical curve with 6=1 and = 4.0 V is the best fit for the experimental verification, 
and we measured an extinction ratio of our apparatus reaching 30 dB. 

In coupling from the IM to the MZI, the actual PBS component has a small rotation 
offset S from the desired 45°, resulting in unequal splitting of the input beam. Thus the real 
device can be modelled as 
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cos 5 — sin 5 cos 5 + sin 6 
- cos 5 — sin S cos S — sin S 



(3) 



In the MZI polarization modulator, the relative phase delay between the two arms is con- 
trolled by high frequency or DC voltages on each modulator of V\ and V2 respectively, in the 
push-and-pull mode. The effect of the PBS and PBC are included implicitly in the Jones 
calculus, as the matrix for the MZI is diagonal: 
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where is the n phase shift voltage for each phase modulator, and $ is the zero-voltage 
phase between two modulators arising from the length imbalance of two arms, denoted AL. 
The value of AL introduces temperature-caused instabilities into the system and plays an 
important role in determining the operating wavelength and voltages of the MZI as seen 
below. To find its value, we perform a wavelength scan over a range of more than 2 nm 
around 1550 nm. When a linear polarizer and a quarter waveplate are placed after the PBC, 
the transmitted intensity I tr exhibits an oscillation: 



hr = (1 + cos (20) cos $ ) • 



(5) 



where $ is the angle of polarizer with respect to the optical axis of MZI, I is the input 
intensity of MZI. $ is a function of wavelength A as $0 (A) = ni^AL. Here rt\ is the 
effective refractive index in the optical fiber, taken as constant over the wavelengths of 
interest; for wider wavelength range, dispersion in the fiber must be taken into account. If 
we use wave number m — jj, then $ becomes a linear function of m as 2ni7rmAL. Eq. 
[5] is plotted in Fig. [3] versus m and compared with the experimental results. Based on the 
data in Fig. [3J we compute AL = 6.0 x 10 -3 m. Since n\ is also a function of temperature 
and wavelength, exactly determining the optimal $ requires further fine tuning the input 
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FIG. 3. The measured intensity oscillation after a linear polarizer and a quarter wave plate when 
the input wavelength of our MZI modulator is varied. The solid line represents simulation data; 
hollow circles are the measured results. 

wavelength. In order to study polarization intensities as Stokes vector elements, we need to 
convert the Jones calculus into a Mueller calculus. The Mueller matrix for a given Jones 
matrix can bd^ converted using M = A - (J ® J*) - A~ l with 

10 l" 
10 0-1 
110 

j -J 

By transforming the matrix [J] in Eq. [2j we find the Mueller matrix corresponding to our 
polarization modulator: 

~1 

cos(6) cos (25) cos(6) sin (25) -sin(0) 
sin(O) cos (26) sin(O) sin (26) cos(6) 
sin (25) -cos (25) 



A 



M, 



o 



(6) 



where 



V1-V2 



7r + $0 is the voltage modulation angle mediated by V\ and V2, which 



are the phase modulator voltages, and $ has been discussed by Eq. ([5]); 5 is the angular 
deviation from perfect 45° input polarization. Now we take the input as linearly polarized 
as in Eq. (1), such that the input Stokes vector is 5 = ^1 1 oj The output of the 



polarization modulator (multiplication by Mq of Eq. rt6fc) is then 
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Eq. indicates the two degrees of freedom that are outside direct control by the phase 
modulator voltages: the intensity imbalance of the arms modelled by deviation angle 5, and 
the initial phase angle $o which stems from the optical path difference between the MZI 
arms. The non-ideal splitting angle S can usually be eliminated by carefully adjusting fiber 
polarization controller FPC_1 in Fig. [2j To set the initial phase $o to the optimal value 
(tt/4) where minimal driving voltages are needed, we varied the wavelength of the telecom 
laser (HP 8168 is illustrated in Fig. [2]) based on Eq. ([5]). After these two adjustments, we 
have the output polarization states 
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Thus, using our modulator, it is possible to output any polarization on the equator of the 
Poincare sphere by controlling the phase voltages V\ and V2. Simulation results are shown 
in Fig. [i] (a). To confirm the accuracy of the simulation, we have performed polarization 
measurements with an Agilent N7788B optical component analyzer; the results of this are 
shown in Fig. [3] (b). To control the phase modulators for these measurements, we supply 
two triangular waves with a n phase difference to the RF ports of the modulators. 

The states of interest for BB84 and corresponding ideal control voltages are listed in 
Table [TJ Driving voltages are minimized when $ = f to reduce power dissipation in the 
driver. The driving voltages are chosen such that the voltage settings are symmetric around 
zero volts, resulting in zero DC bias voltage. The maximum sampling speed of the polar- 
ization analyzer we used is 0.1 MHz, much less than the 76 MHz system clock rate. Thus, 
to characterize the polarization states at higher speed, we developed a new method which 
we describe below. 
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TABLE I. The driving voltage settings of two arms and corresponding output polarization states. 
The value of voltages labeled in the right hand figure is V\ —V2. 
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FIG. 4. Poincare sphere traces of the light passing through the polarization modulator, due to 
a triangular shaped voltage, (a) Simulation results; (b) experimental measurements, (a) and (b) 
are determined after the polarization controller FPC_2 and include the unitary transformation 
QWP(|) to bring the modulation to the equator of the Poincare sphere. The polarization states 
are measured by an Agilent N7788B Optical Component Analyzer, and the Polarization Navigator 
software package. 



D. Polarization state characterization 

We wish to know the exact output polarization state described by Eq. ^ at full mod- 
ulation speed. It is impossible to do a direct characterization of the Stokes vectors using 
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standard polarimeters due to the high repetition rate of the modulator. Here we use a quar- 
ter waveplate (QWP) followed by a linear polarizer with a repeated modulation sequence 
that allows us to extract the Stokes parameters of each state in the sequence. The intensity 
Sf after the QWP and polarizer is calculated in RefP^ 

Sf =\{S + (Si cos 2 fa + S 2 sin 2 fa) cos 2 (a 3 - fa) 
+ [(£2 cos 2/3 j — Si sin 2 fa) cos A + S3 sin A] 

x sm2(a j -fa) }, (9) 

where aj is the angle of the linear polarizer's transmission axis and fa is the angle of the fast 
axis of the QWP — in our experiment, angles aj and fa are both with respect to horizontal. 
The retarder phase delay A is tt/2 for a QWP, however in real devices there exists non-trivial 
deviations to this ideal value. In our setup, we use a QWP from Thorlabs and find 7% offset 
from |. 

At least three measurements with the settings defined in Table [Til were performed to obtain 
three components (S^, S^, and S3) taken at 2 MHz modulation speed. The results are 
displayed in Fig. [5] (a). This characterization scheme utilized an optical detector (Thorlabs 
PDA10CF) which had a bandwidth of 150 MHz. 

In Fig. [5] (a), there is a time displacement of around 100 ns between the differential voltage 
V1 — V2 and polarization states due to the delays between the optical and electrical signals. 
The spike overshoot is due to the finite transition time across multiple polarization states, 
when the differential driving voltage has the largest step. The transition ring ripples after 
the spike is caused by the limited response time of the optical PIN diode, but the actual 
polarization states should not be affected. Moreover, since the upconversion SFG process 
only occurs in a time overlap (less than 1 ns) when both telecom and 810 nm photons arrive, 
the polarization states of telecom photons within a few percentage of each repetition cycle 
are responsible for the green photon states. We are able to completely avoid the transition 
spike by moving the overlap time to the end of each cycle. 

In the subsequent simulation, we have modeled the entire IPM system as well as the po- 
larization characterization optics. The imperfection of QWP was also take this into account. 
The degree of polarization (DOP) was experimentally found to be over 99%, hence So can be 
determined by the averaged maximum intensity. We then use the relations Si = 2S± — So, 
S2 = 2S2 — So, and S3 = 2S3 — So to determine the exact polarization states in H/V, D/A, 
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TABLE II. Measurement settings for polarization state characterization 
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and L/R bases, respectively, for each modulation setting (Vi and V2). The simulation results 
compared with the extracted Stokes vector elements Si, £2 and £3 (in hollow square) are 
displayed in Fig. [5] (b). For each extraction point, the x error bar indicates the length of 
time (125 ns) to average the polarization states; the y error bar indicates the uncertainty of 
each extraction. 

E. Temperature Stability of the Modulator 

The temperature stability of the polarization state plays a critical role for long-term 
operation of this device. To measure the performance of a given polarization state against 
the variation of ambient temperature, we use a constant optical power and wavelength at the 
input port of our polarization controller, and measure the output power of four polarization 
components by our in-house-developed polarimetry method described in the previous section. 

Fig. [6] shows the measurement results for H/V and D/A components. The temperature 
is maintained by a Thorlabs TC-200 controller. In this measurement, no voltage is applied 
to either the phase or intensity modulators, but the initial phase $ in the MZI becomes a 
function of temperature. By carefully eliminating the L/R components, the output polariza- 
tion state sweeps across the equator of the Poincare sphere (provided that the temperature 
variation is big enough) as indicated by Fig. |6| 

However, at higher temperature settings when the temperature setting is > 30°C, the 
temperature actuator (heater) cycles through on and off states more frequently, resulting 
in strong fluctuations and poor visibility. This is due to the simple point contact heater 
used — it takes some time for the temperature of the arms of the MZI to stabilize, and 
more importantly two arms are not heated up equally leading to adversely increasing the 
imbalance of MZI. A two-stage heater would be desirable for achieving better long term 
temperature stabilization. 
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FIG. 5. (a) The oscilloscope traces to measure Stokes vector including differential voltage V\ — V% 
between two phase modulators located at two arms of MZI; three measurements for 5+, S£ and S% 
at a wavelength of 1550.5 nm. All of the raw signal traces were collected by an Agilent DSO8104A 
Innniium Oscilloscope. Each horizontal division was set to be 500 ns in the plot; the differential 
voltage plot had 2 volts per division for its vertical trace; while the rest three plots had 100 mV 
per division, (b) The simulation results for wavelength A = 1550.5 nm, with measured three Stokes 
vector elements normalized to So (in blue hollow squares with error bars). The simulation includes 
the imperfect retardance (7%) of the QWP; the time delay between differential voltage to the 
polarization states has not been taken into account. 



III. EXPERIMENTAL RESULTS OF INTENSITY AND POLARIZATION 
MODULATOR IMPLEMENTING BB84 PROTOCOL WITH DECOY 
STATES 

A. Theory of BB84 protocol with decoy states 

WCP QKD, attractive for its simple designs and high rates, is subject to the photon 
number splitting attackP^ due to the nonzero probability of producing multiple photons in a 
single laser pulse. To securely achieve any reasonable transmission distance requires the in- 
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FIG. 6. The measured power of four polarization states as temperature varies. The setting of 
temperature is monitored and maintained by a Thorlabs TC-200 temperature controller driving a 
point contact heater in the enclosure of the MZI polarization controller. 

elusion of decoy states, i.e. pulses with differing intensities used to bound the eavesdropper's 
information from multi-photon eventd^Sl. We implement a polarization-encoded protocol 
with vacuum+weak decoy stated. To calculate the final secure key rate, we must estimate 
the single-photon gain and error rate. We calculate a lower bound of single photon gain Q\ 
from Ref p2 as 

Q\ = (q v <? - Q,^- 2 - ^^Y ) , (10) 

where ji is the average photon number for signal states and v for decoy states, and Q^/v 
are the gains for these two states. 1q, the vacuum yield, is determined by the detector dark 
counts and background noise, and was measured between non- vacuum pulses. Similarly, the 
upper bound on the error rate, , of single photon states is 
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where E y is the total error rate for decoy states, and e = 0.5 is the vacuum error rate. 



After evaluation of the parameters in Eqs. (10) and (11), we can find the lower bound of the 
secret key rate per pulse, R: 



R = qL Kul) {-Q,f (QBER) H 2 (QBER) + Q\ [l - H 2 ( £ f )] } . 



(12) 



where q — 1/2 is the basis reconciliation factor, = N /jL /(N IJL + N u ) is the ratio of signal 
pulses to the total pulses, and QBER is the total quantum bit error rate for signal pulses. For 
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FIG. 7. (a) Experimental measurements for the decoy+signal QKD states accumulated in the 
histogram format (se^Sl in which Bob's detector was fixed to an H/V basis measurement.). 1 45) 
indicates either a D or A polarization state which is not resolved by our two-detector receiver in 
this basis, (b) Stability measurement of the overall QKD system-wide QBER over consecutive 160 
minutes at 25 dB total loss, showing the averaged QBER at 1.8 ± 0.9%. The system clock rate is 
76 MHz triggered from the Ti-Sapphire mode-locked laser. 

the calculations here, a constant error correction efficiency /(QBER) of 1.22 was assumed. 

B. Experimental Demonstration of High Loss QKD System 

Our complete system exhibits high performance: when the 810 nm laser is not mode- 
locked (i.e. continuous wave), the output photons at 532 nm have a fidelity >99% with 
the desired polarization states. Upon mode-locking and fast modulation, we find that our 
system fidelity including all losses and error sources is maintained at >98% in the H/V 
polarization basis, and >95% in the D/A basis. The modulation achieves accurate average 
photon number per pulse for both decoy and signal states, sufficient to perform successful 
QKD. Fig. [7] (a) demonstrates the measured photon output for signal and decoy states with 
various polarizations. We additionally implemented an automated polarization alignment 
procedure in the quantum receiver to keep high fidelity for all four states, the details of 
which will be published elsewhere. Fig. [7] (b) illustrates a QBER stability measurement over 
more than 160 minutes with ambient temperature fluctuation <0.5°C. The QBER shows 
an average value of 1.8 ± 0.9%, a very small variation. Our simulations of ground-to-space 
quantum channels show that the overall QBER of our system is sufficient for a satellite 
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FIG. 8. The projected performance of our decoy state BB84 QKD system. We show the simulation 
results of QBER, secure key rate probability per pumping laser pulse, and the final key rate under 
total channel loss. We assume two pumping laser rates: the currently used Ti-Sapphire laser 
pumping rate of 76 MHz (solid red line), and the projected 1 GHz clock rate (dashed blue line). 



uplinkP. 

Finally, we calculate the secure key rate based on Eq. (12). The result is shown in Fig. [8j 
in which the detector efficiency is 60% and the dark count rate is 50 counts/s. The intrinsic 
QBER is 1%, which accounts for the optical misalignments in the telecom modulator, SFG 
source, and receiver. In practice, stray light will lead to slightly higher measured QBER 
than the simulated value under the same channel loss. The signal and decoy states are 0.6 
and 0.2 photons per pulse respectively. QKD beyond 60 dB total loss is possible with this 
system. 
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IV. CONCLUSION AND OUTLOOK 



We have presented our design and implementation of a Mach-Zehnder configuration, all- 
fiber, intensity and polarization modulator. This design features very high switching speed 
and high polarization visibility. Our modulator serves as a secret key encoder which can 
easily be upgraded to GHz speeds owing to the rich choice of photonic devices for telecom 
wavelengths. We used this modulator and the SFG process to produce 532 nm green photons 
at a repetition rate near one hundred MHz, which is the highest speed for this operating 
wavelength demonstrated so far. The green photons show excellent polarization state fidelity, 
and are compatible with the highest figure of merit single photon detectors commercially 
available. 

The quantum optical simulations in Fig. [8] show QKD is possible over channels beyond 
60 dB of loss. Additionally, we measured the raw key rates, average photon numbers for 
signal and decoy states, and quantum bit error rates with high channel losses. The perfor- 
mance fulfills the design goals of the telecom modulator and agrees with our quantum optics 
simulation results. Assuming an average loss of 45 dB for a LEO satellite passage, the final 
secure key rate is 155 bit/s. 

If we employ a 1 GHz repetition rate pumping laser, the operating clock rate of our tele- 
com modulator can be readily boosted accordingly, thanks to the technological advancement 
of FPGA^ and RF power electronics^. Assuming this increased speed, the QBER, secure 
key rate per pulse and per second are displayed in Fig. [8] as dashed lines. At the same 
average LEO satellite channel loss, the secure key rate will be ^2300 bit/s, an increase of 
more than an order of magnitude, and better than the repetition rate increase factor owing 
to an improved signal-to-noise ratio. Consequently, we present our system that is able to 
overcome the high channel loss that plagues proposals for a satellite uplink QKD system 
with our current setup. Our experimental results verify the feasibility of using such a source 
in a near term satellite mission. Using higher speed modulation such as 1 GHz rate, it is 
possible to achieve higher secrete key rate and to withstand higher channel losses. 
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